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ABSTRACT 

Recent spectroscopic observations of galaxies in the Fornax-Cluster reveal nearly un- 
resolved 'star-like' objects with red-shifts appropriate to the Fornax-Cluster. These 
objects have intrinsic sizes of « 100 pc and absolute B-band magnitudes in the 
range — 14 < Mb < —11.5 mag and lower limits for the central surface brightness 
^ 23 mag/arcsec^ (Phillipps et al. 2001), and so appear to constitute a new popu- 
lation of ultra-compact dwarf galaxies (UCDs) . Such compact dwarfs were predicted 
to form from the amalgamation of stellar super-clusters (Kroupa 1998), which are rich 
aggregates of young massive star clusters (YMCs) that can form in collisions between 
gas-rich galaxies. Here we present the evolution of super-clusters in a tidal field. The 
YMCs merge on a few super-cluster crossing times. Super-clusters that are initially as 
concentrated and massive as knot S in the interacting Antennae galaxies (Whitmore 
et al. 1999) evolve to merger objects that are long-lived and show properties compa- 
rable to the newly discovered UCDs. Less massive super-clusters resembling knot 430 
in the Antennae may evolve to uj Cen-type systems. Low-concentration super-clusters 
are disrupted by the tidal field, dispersing their surviving star clusters while the re- 
maining merger objects rapidly evolve into the /iB — Mb region populated by low-mass 
Milky- Way dSph satellites. 

Key words: galaxies: interactions - galaxies: formation - galaxies: star clusters - 
galaxies: dwarfs - globular clusters: u: Cen (NGC 5139) - methods: N-body simulations 



1 INTRODUCTION 

In a recent spectroscopic survey of all objects around 
NGC 1399 in the Fornax-Cluster, carried out to find all 
dwarf galaxies, Phillipps et al. (2001) report five ultra- 
compact dwarf galaxies (UCDs; see also Hilker et al. 1999). 
These have spectra typical for late-type metal-rich and 
old stellar populations, are marginally resolved and have 
red-shifts comparable to the Fornax-Cluster, ruling out ei- 
ther faint background galaxies or foreground stars. The 
UCDs have intrinsic sizes of around 100 pc, absolute B- 
band magnitudes of —14 < Mb < —11.5 mag and lower 
limits (due to failure to resolve their cores) of the cen- 
tral surface-brightness ~ 23 mag/ arcsec^ . Further anal- 
ysis of photographic plates show no sign of low-luminosity 
envelopes around these objects which rules out the possi- 
bility that they are nucleated dwarf ellipticals (dE,N) with 
faint envelopes. Therefore, these objects are either extremely 
compact dwarf galaxies or extremely large and massive 
(10^ — 10* Mq, assuming M/Lb = 3) star clusters. In the 
^B — Mb diagram they fall in the empty region between 
'ordinary' dwarf galaxies and globular clusters (Fig. 

High-resolution HST-images of the star forming regions 
in interacting galaxies like the Antennae (Whitmore et al. 
1999; Zhang & Fall 1999) or Stephan's Quintet (Gallagher 



et al. 2001) resolve some of these regions into dozens to hun- 
dreds of young massive star clusters. According to Whitmore 
et al. the individual clusters have efi'ective radii of about 
4 pc and masses lO^-lO® M©, with a mass-spectrum follow- 
ing a steep power law *a/ oc M"^ (Zhang & Fall 1999). The 
striking point is that these young star clusters are them- 
selves clustered into groups of a few to several hundred star 
clusters spanning projected regions of a few 100 pc, with a 
higher cluster concentration at the centre. Measurements of 
the relative velocities between the star clusters within such 
objects are now becoming available. Preliminary results in- 
dicate ~ 20 km/s (B. Whitmore, private communication), 
which is consistent with virial masses « 10^ Mq . Age deter- 
minations in the Antennae show that these star clusters are 
extremely young (3-7 Myr). While in the Antennae young 
massive star clusters are preferably found in the central re- 
gion of the interacting pair, NGC 7319 in Stephan's Quintet 
has young star clusters in the long tidal arm and the intra- 
group region north of NGC 7319 (Gallagher et al. 2001). 

Examples of clusters of star clusters are shown in Fig. ^ 
as thick crosses. Of these, knot "S" is particularly remark- 
able, as its mass is estimated to be about 2.1 x 10^ 
assuming it has a mass-to-light ratio in the photometric V- 
band of M/Lv = 0.076 for an age of 10 Myr (fig. 7 in Smith 
& Gallagher 2001) and Mv = -15.8 (section 4.4.2 in Whit- 
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Figure 1. Central surface brightness, ^o,Bi vs absolute photometric B-band magnitude, Mb (the Kormendy diagram). The solid circles 
are data for Milky-Way globular clusters (GC, Harris 1996), Local-Group dwarf galaxies (dirr, Andl and dSph, Mateo 1998), elliptical 
galaxies (E, Peletier et al. 1990). The positions of disk galaxies are indicated by "Spirals" (cf. Ferguson &: Binggeli 1994). The newly 
discovered UCDs are shown by arrows with lower-limits on /ib (Phillipps et al. 2001). Three "knots", (KnS, Kn430, Kn225, tables 1 
and 2 in Whitmore et al. 1999), observed in the interacting Antennae galaxies are shown as thick crosses. Knots KnS and Kn430 are 
roughly 10 Myr old so that M/Lv ~ 0.076 (fig. 7 in Smith & Gallagher 2001), while Kn225 is about 500 Myr old. The fading of all 
three knots to M/Lg = 3 due to stellar evolution but assuming no morphological evolution (eq. is indicated by the dashed arrows. 
The thick six-pointed star shows the merged super-clusters predicted by Kroupa (1998) to result from the dynamical evolution of such 
knots, which consist of dozens to hundreds of young star clusters. 



more et al. 1999). Two other knots (Kn430 and Kn225) are 
plotted in Fig. hi with fading vectors to M/Lb = 3 to in- 
dicate their position after a few Gyr if the shape of their 
density profile remains unchanged, 

~ constant, (1) 

where and Ltot are the central and integrated luminosity, 
respectively. 

The dynamical evolution of these knots is an interest- 
ing problem to consider because the merging of many star 
clusters is not well understood but is a pre-requisite for the 



understanding of the future fate of such knots. By scaling 
the numerical results of Garijo, Athanassoula & Garcia- 
Gomez (1997, hereinafter GAG), who studied the forma- 
tion of cD galaxies from clusters of 50 equal-mass galax- 
ies, to the present problem of a super-cluster containing 50 
star clusters, each cluster having a mass Mci = 10^ Mq, 
Kroupa (1998) studied the dynamical evolution of such "stel- 
lar super-clusters". These are not to be confused with the 
commonly-used terminology 'super-star-cluster', SSC, which 
is one young star cluster with a mass typical for globu- 
lar clusters (e.g. Smith & Gallagher 2001). Depending on 
the initial concentration of such an object, the majority of 



Formation of UCDs 3 



star-clusters merge within a few super-cluster crossing times. 
Within 100 Myr a stellar system forms which has a relax- 
ation time significantly longer than a Hubble time. It is thus 
a low-mass spheroidal dwarf galaxy (as opposed to the well- 
known dwarf-spheroidal satellite galaxies, denoted dSph), 
with a half-mass radius of between about 45 to 100 pc, pos- 
sibly with a high specific frequency of globular clusters if 
some survive the merging phase. From figs. 16-18 in GAG 
the central surface density for essentially all models (col- 
lapsing super-clusters and super-clusters initially in virial 
equilibrium) is 1 mass-unit/length-unit^. This translates to 
/i — lOOOMo/pc^ and to the apparent-magnitude surface 
brightness /ib ~ 19.48 mag/arcsec^ using 



/iB[mag/arcsec ] = 



-2.51ogio (m^) 

Mb,0 -5 + 5, 



(2) 



where Mb,© = -1-5.41 is the absolute B-band magnitude of 
the Sun, and S = 26.572 = 2.51ogio/^ with / = 2.06 x 
10^ arcsec per radian, and assuming the B-band mass-to- 
light ratio M/Lb = 3 in solar units. Assuming that 50- 
100 per cent of all clusters merge to form the spheroidal 
dwarf galaxy, its integrated absolute B-band magnitude, 

AfB = - 2.51ogio (m^) +Mb,0, (3) 

where M is the mass of the object in Mq, lies in the range 
-12.64 < Mb < -11.89. 

FigTQ shows that the predicted spheroidal dwarf galax- 
ies fall in the region of the Kormendy diagram that the UCD 
data of Phillipps et al. (2001) cover, nicely verifying the 
prediction. The hypothesis that the UCDs were formed as 
super-clusters appears to be a reasonable possibility, since 
gas-rich galaxies will have merged in profusion in sub-groups 
during the early phases of the formation of the Fornax 
galaxy cluster. The UCDs may thus have formed during 
the coalescence of a subgroup of late-type gas-rich galaxies, 
which was infalling into the Fornax cluster. Tidal shocking of 
such UCDs and/or initially less-massive super-clusters may 
account for some if not the majority of known spheroidal 
dwarf galaxies found in large numbers in rich galaxy clus- 
ters (Gallagher, Conselice & Wyse 2001). 

Kroupa (1998) studied isolated super-clusters, but tidal 
shaping of the merger objects may be significant. In this 
contribution we report self-consistent stellar-dynamical cal- 
culations of the merging process in super-clusters whilst they 
are immersed in the tidal field of a parent galaxy. A detailed 
study of the merging time-scales and efiiciencies for super- 
clusters on circular orbits shows that the star clusters merge 
very efficiently on a few dynamical time-scales of the super- 
cluster (Fellhauer et al. 2001). In this paper we focus on a 
comparison with the recently discovered UCDs and perform 
computations of super-clusters with varying initial concen- 
trations on eccentric orbits. 

The models are explained in Sect. |^ where a brief de- 
scription of the N-Body code is given. In Sect. |^ the prop- 
erties of the merger objects on circular orbits are brieffy 
described and a model on an eccentric orbit is discussed in 
detail. The evolution of a low-mass extended merger object 
is also studied. Conclusions are presented in Section 0. 



2 NUMERICAL METHOD AND SETUP 
2.1 The code 

The orbital integration of the particles is performed with 
the particle- mesh code Superbox (Fellhauer et al. 2000). 
In Superbox densities are derived on Cartesian grids us- 
ing the nearest-grid-point (NGP) scheme. From these den- 
sity arrays the potential is calculated via the fast Fourier- 
transformation. Forces are obtained using higher-order dif- 
ferentiation based on the NGP scheme but comparable in 
precision with standard CIC (cloud-in-cell) algorithms. The 
particles are integrated using a fixed time-step Leap-Frog 
algorithm. For an improved resolution at the regions of in- 
terest, a hierarchical grid-architecture with two levels of high 
resolution sub-grids are used for each star cluster in the 
super-cluster. These sub-grids track the density maxima of 
the individual star clusters, and are adjusted at the begin- 
ning of the computation to meet individual requirements. 

The sizes of the two high resolution grid-levels are cho- 
sen such that the innermost grids with the highest resolution 
(1 pc per grid-cell length) cover the star clusters, while the 
medium-resolution grids cover the merger object to ensure 
that the forces between each cluster are treated at least at 
medium resolution. Finally, the outermost grid, which stays 
fixed and has the lowest resolution, covers the whole orbit of 
the super-cluster around the parent galax;y. The time-step 
used in the calculations is dt = 0.1 Myr, which is signif- 
icantly shorter than the crossing time in the star-clusters 
{tcr « 6 Myr). 

Out notation is as follows: we use large letters for galac- 
tocentric coordinates {X, Y, Z, R) and derived quantities, 
and small letters {x,y,z,r) for super-cluster-centric coor- 
dinates and derived quantities. 



2.2 Orbits 

The super-cluster orbits within an external potential of a 
parent galaxy given by 

1, 



HR) 



-V^r-ln(i?^,i + i?^ 



(4) 



with core radius _Rgai = 4 kpc and circular velocity Vo = 
220 km/s. Dynamical friction is also analytically added using 
the Chandrasekhar expression (Binney & Tremaine, their 
eq. 7-18): 



AVm 
At 



47rln(A)G^£igaiMs, 



(5) 



erf 



V2Vm 



■ exp 



Vm 



Vm 



with ln(A) = 1.5 (Velazquez & White 1999; Fellhauer et al. 
2000) the so-called Coulomb-logarithm, q^^i is the density of 
the galactic halo at the position of the super-cluster, M^c is 
the total mass of the super-cluster, Vm is the orbital velocity 
of the super-cluster in the tidal field of the host galaxy and a 
is the velocity dispersion in the galactic halo at the position 
of the super-cluster. 

We adopt the tidal field of a Milky- Way type galaxy, 
but our results are applicable to a galaxy cluster such as For- 
nax in the sense that super-clusters that remain stable over 
many highly eccentric orbits will also survive in the potential 
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Table 1. Properties and orbits of the model super-clusters. The columns give the following data: name of the 
model, number of star clusters in the super-cluster, Plummer-radius of the super-cluster, crossing time of the super- 
cluster, total mass of the super-cluster, perigalacticon distance, apo-galacticon distance, orbital period, concentration 
parameter, strength of tidal field parameter. 



name 




^sc 

pi 


J-SC 

''cr 




^peri 


-Rapo 


^orb 


a 








[pc] 


[Myr] 


[IO^Mq] 


[kpc] 


[kpc] 


[Myr] 






SD05A 


20 


50 


7.4 


2.0 


10 


10 


300.7 


0.12 


0.62 


SD06A 


20 


75 


13.5 


2.0 


10 


10 


300.7 


0.08 


0.93 


SD09A 


20 


50 


7.4 


2.0 


20 


20 


569.6 


0.12 


0.42 


SDIOA 


20 


75 


13.5 


2.0 


20 


20 


569.6 


0.08 


0.60 


SDll 


20 


150 


38.3 


2.0 


20 


20 


569.6 


0.04 


1.28 


SD17B 


20 


50 


7.4 


2.0 


50 


50 


1400.7 


0.12 


0.23 


SD18A 


20 


75 


13.5 


2.0 


50 


50 


1400.7 


0.08 


0.35 


SD19E 


20 


150 


38.3 


2.0 


50 


50 


1400.7 


0.04 


0.70 


R06 


20 


300 


108.4 


2.0 


30 


60 


900 


0.02 


2.0 


EOl 


32 


50 


9.4 


1.2 


2 


20 


250 


0.08 


1.29 



well of Fornax, whereas super-clusters that are torn apart 
immediately are not likely to survive in a galaxy cluster. 
Gallagher, Conselice & Wyse (2001) discuss dwarf-galaxy 
survival in more detail, and we note that spheroidal dwarfs 
without dark matter can readily survive passages through 
the centres of galaxy-clusters owing to the large encounter 
velocities. We do not aim to model the UCDs observed in 
the Fornax cluster in detail, which is impossible in any case 
given that little is known about their orbits, but merely 
concentrate on the generic properties of the merger objects 
when they are subject tidal fields that vary in strength. 

To study the merging time-scales and efficiencies circu- 
lar orbits are adopted where the centre of the super-cluster 
is placed at a distance R — D from the galactic centre. The 
distance D is varied in the range of 5-100 kpc. A detailed 
discussion of the involved time-scales can be found in Fell- 
hauer et al. (2001). To study the evolution of the merger 
objects in a time- dependent tidal field the parameter space 
is extended to eccentric orbits. 



2.3 Super-clusters 

A range of models (Table 0) are setup to investigate the 
importance of the tidal field and initial concentration of the 
super-cluster. 

The young massive star clusters are modelled as 
Plummer-spheres with Plummer-radii Vpi = 4-6 pc and cut- 
off radii 20-30 pc. The clusters all have masses of 10® Mq 
if no mass spectrum is taken into account, or 10® Mq and 
3.2 • 10^ M0 to model a mass spectrum. Estimating the re- 
lax;ation time frciax of these star clusters (Spitzer & Hart 
1971), 



= 0.138- 



(6) 



\/Gln(0.4 n) 

where M^i is the mass of the star cluster, is the half- 
mass radius, the average stellar mass (m) — 0.4 M0 (using 
the universal IMF, Kroupa 2001) and n the number of stars. 
This leads to troiax ~4 Gyr for the more massive star clusters 
and 2.5 Gyr for the less massive ones. These time-scales 



are much longer than the cluster-cluster merging time-scale. 
The resulting merger objects have relaxation times longer 
than a Hubble time. Surviving star clusters are only counted 
and their position is traced; their internal evolution is not 
studied here. Therefore it is possible to carry out the orbital 
integration with a particle-mesh code which suppresses two- 
body relaxation processes. 

To model a super-cluster, A'o, of these clusters are 
placed in a Plummer distribution with Plummer radius (r'^i) 
ranging from 50-300 pc and cut-off radius (r^Jjt). The super- 
cluster is initially in virial equilibrium. 

To describe our parameters with dimensionless variables 
we define 



/3 = 



/ sc 

= rpi/rpi, 

lit /'"tidal J 



(7) 
(8) 



where rtidai is the tidal radius of the super-cluster. In the 
case of eccentric orbits the tidal radius is estimated at peri- 
galacticon, 

1/3 



nidal = 



3Mgal(i?) 



Rr, 



(9) 



While a describes how densely the super-cluster is filled with 
star clusters, /3 indicates the strength of the tidal forces act- 
ing on the super-cluster. With our choices of the parameters 

and D, a ranges from 0.02 to 0.12 and /3 falls in the 
interval 0.2 to 2.7. 

The relaxation time of the cluster system is very short, 
i.e. using Eq. ^ with n — No would lead to a nominally 
shorter relaxation time than the crossing time of the super- 
clusters (about 0.3-0.4 tl'^). This shows that such systems 
are heavily collision-dominated and therefore the merging 
process is violent and the merging timescales should be very 
short. 



3 RESULTS 

As discussed by Kroupa (1998) and Fellhauer et al. (2001) 
the evolution of the young super-cluster is dominated ini- 
tially, during the violent merging time, by cluster-cluster 
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Table 2. Sizes of the merger objects, measured at t - 
the bound mass is taken into account (cf. Fig. ^). 



1 Gyr and at t = 5 Gyr for different choices of a and /3. Only 



name 


a 


/3 


Half-Mass 


Radius [pc] 


90 %-Mass Radius [pc] 


Tidal Radius [pc] 








t = 1.0 Gyr 


t = 5.0 Gyr 


t = 1.0 Gyr 


t = 5.0 Gyr 


t = 1.0 Gyr 


t = 5.0 Gyr 


SD05A 


0.12 


0.62 


39.0 


40.5 


123.2 


135.4 


400 


400 


SD06A 


0.12 


0.42 


48.0 


50.2 


151.5 


150.9 


600 


600 


SD09A 


0.12 


0.23 


38.3 


39.0 


125.5 


124.2 


1,120 


1,120 


SDIOA 


0.08 


0.93 


79.2 


77.3 


200.0 


188.9 


390 


390 


SDll 


0.08 


0.60 


43.0 


44.5 


146.1 


150.9 


610 


610 


SD17B 


0.08 


0.35 


71.2 


72.4 


230.4 


221.5 


1,110 


1,110 


SD18A 


0.04 


1.28 


181.2 


144.6 


491.1 


348.1 


580 


580 


SD19E 


0.04 


0.70 


236.6 


163.2 


1,058.8 


576.5 


1,010 


910 


R06 


0.02 


2.00 














EOl 


0.08 


1.29 


54.6 


41.9 


149.2 


146.3 


188 


152 



merging and interaction processes. The merging efficiency is, 
in general, very high. As long as the cutoflt radius is smaller 
than the tidal radius (/3 < 1.0) almost all star clusters merge. 
The total number of clusters that merge decreases linearly 
and slowly with increasing f3 > 1. 

3.1 Circular orbits 

The merger objects consisting initially of A'o ~ 20 massive 
star clusters (each with 10^ M©) on circular orbits have half- 
mass radii and sizes summarised in Table ^ The time-scale 
for the merging process is very short. Within 5 — 10 x tl'j. 
{tl'j. being the internal crossing time of the super-cluster) the 
initial violent merging phase is over. During this phase the 
number of remaining clusters decreases exponentially, the 
decay thus being similar to the radioactive decay law. For 
initially concentrated super-clusters {a > 0.08) this phase 
lasts ;$ 100 Myr. The resulting merger object is very dense 
and stable over 10 Gyr. In Fig. |l| it coincides with the Kr98 
prediction. 

As long as 0.08^ a 0.12 there is almost no evolution 
in mass and size of the merger object. If a 0.04 or /3 > 1.0 
the merger object looses a part of its mass and shrinks until 
it reaches a stable configuration, and is associated with a 
population of surviving star clusters that spread out along 
the orbit. The merger object lies close to the Kr98 prediction 
in Fig. ^ despite the mass lost. 

Between 70 and 95 % of the mass of the super-cluster 
ends up in the merger objects, which have central densi- 
ties of about 650 Mq/pc^ and exponential scale lengths 
of about 9 pc. The surface density can be well fitted by 
a de-Vaucouleur-profile, the line-of-sight velocity dispersion 
is about 20 km/s. The shape of the objects is spherical with 
axes-ratios larger than 0.9. 

3.2 Dense super-cluster on an eccentric orbit 

The models are extended to include eccentric orbits to study 
the infiuence of tidal heating acting on the super-clusters 
and the resulting merger objects (Table |^). As reference 
model in this paper we start with a dense super-cluster con- 
taining A^o = 32 star clusters, 3 having a mass of 10*^ Mo 



and 29 having 3.2- 10^ M© to mimic a mass spectrum of the 
form oc . The star clusters have a Plummer radius 
of rpi = 4 pc and crossing times of 0.75 and 1.32 Myr, respec- 
tively. The super-cluster they are distributed in, has a total 
mass Msc = 1.23 x lO'^ Mq, a Plummer radius = 50 pc, 
a cut-off radius of r^ut = 250 pc (giving a = 0.08) and 
a crossing time of 9.4 Myr. The super-cluster is placed at 
apo-galacticon (D — 20 kpc) on a highly eccentric orbit 
with eccentricity e — 0.82, where 

e = ^apo - Rpcri ^^g^ 

At perigalacticon /Jpori = 13 = 1.29. However, since the orbit 
starts at apo-galacticon and has an orbital period Torb ~ 
250 Myr the super-cluster has /3 < 1 for a time-span nt x f^r, 
nt ~ 10, giving enough time for the super-cluster to evolve 
essentially as the models studied above. Most star clusters 
merge within 5 x t^^ (Fig. The number of star clusters 
decreases exponentially, 

N{t) = iVo.exp(-^j, (11) 

with t measured in units of the crossing time of the super- 
cluster and T = 5 tcr which corresponds to 47 Myr for model 
EOl. This exponential decrease is shown as a straight line 
in Fig. 1^ The initial violent merging phase is dominated by 
clusters on highly eccentric orbits within the super-cluster. 
Star clusters initially on near-circular orbits in the super- 
cluster survive until ^ 5 tcr , but ultimately they merge due 
to dynamical friction within the super-cluster. 

The mass of the merger object is Mmo = 1.12 • lO'^ M© 
which corresponds to 91.5% of the total mass of the super- 
cluster. The remaining mass is lost due to internal heating 
from cluster-cluster collisions and subsequent loss due to 
the tide. The mass of the merger object decreases with every 
perigalacticon passage. The averaged mass-loss rate is 5.7 • 
10^ Mo/Gyr. After 10 Gyr the total mass of the remaining 
object is 0.59 • lO'' M©. 

By t = 1 Gyr, the gravitationally bound particles com- 
prising the merger object have a diameter of several 100 pc 
with an elliptical envelope (e ~ 0.6) but a spherical and very 
dense core (Fig. ^). Bound particles have negative binding 
energies relative to the merger object. The half-mass radius 
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Figure 2. The number of remaining star clusters, A^rcm, vs. time 
for model EOl. Within the first 5 crossing times of the super- 
cluster the merging can be fitted with an exponential with an 
exponential decay time of about 47 Myr. This is the same fit as 
already obtained by Fellhauer et al. (2001) for a large variety of 
super-clusters. 



eOl E(, = 16.7 mag/arcsec»»2 20 pc Resol. I = 600 Myr 

; ' ^ 




-3.3 -3 -2.8 -2.6 -2.4 -2.2 

X [kpc] 

Figure 3. Contour plot of the merger object EOl at 600 Myr 
shortly after the violent merging phase is over. While the central 
area is circular, the envelope shows an ellipticity of about 0.6. 




L [Gyr] 

Figure 4. Lagrangian radii containing 10, 20, 90 % of the 
mass of merger object EOl vs. time. Solid lines are all particles, 
dashed lines are the bound particles only. 




1 10 100 

r [pc] 



Figure 5. Radial density profile of merger object EOl at t = 1 
(crosses), 5 (points) and 10 Gyr (circles). Profiles can be fitted 
by two exponentials (dashed lines, according to eq. |l2[ ) with the 
overall fitting function as solid line. Vertical lines show the tidal 
radii at t = 1 and 10 Gyr derived from Eq. M. 



rh = 54.6 pc and r{90%) = 149.2 pc holds 90 % of the mass 
(Fig. Examining the surface density profile (Fig. ^), the 
core radius rcorc = 10 pc (the radius where the surface den- 
sity drops to half the central value). The tidal radius at 
perigalacticon, estimated from Eq. ^ is rtidai ~ 188 pc. The 
morphology of the object stays almost constant during the 



whole computation time. At t = 10 Gyr the half-mass radius 
is rh = r(50%) = 41.9 pc, r(90%) = 146.3 pc, r.orc = 9.5 pc 
and rtidai = 152 pc, as shown in Fig. ^. 

The structure of the merger object (MO) can be best 
fitted by a double exponential profile. The MO shows a very 
dense core surrounded by an exponential envelope. Defining 
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Figure 6. Surface density profile of tlie merger object EOl. Pro- 
files are displayed at t = 1 (crosses), 5 (dots), 10 Gyr (cir- 
cles). Solid thick line shows a de-Vaucouleurs fit to the data at 
t = 10 Gyr. Also shown as dashed lines are two exponentials as a 
fitting function with the sum of both shown as thin solid curve. 
Right vertical dotted line shows the tidal radius taken from Eq. ^ 
left vertical dotted line shows the core radius (both at t = 10 Gyr). 
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Figure 8. Time evolution of the central values of space density, 
surface density and line-of-sight velocity dispersion of merger ob- 
ject EOl. 
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Figure 7. Line-of-sight velocity dispersion profiles of merger ob- 
ject EOl at t = 1 (crosses) and 10 Gyr (boxes). Profiles can be fit- 
ted by exponentials. Fitting parameters for t = 1 Gyr: inner: (tq = 
12.7±0.1 km/s, rexp = 110± 10 pc; outer: ctq = 12.4±0.2 km/s, 
r^xp = 222 ± 5 pc; t = 10 Gyr: inner ctq = 11.6 ± 0.6 km/s. 



: 90±2 pc; outer: 



10.5±0.3 km/s. 



185 ±10 pc. 



the profile as 
p(r) = po exp 



(12) 



the central density at t = 1 Gyr is go = 1030 ± 40 M© /pc^ 
and the exponential scale length of the core is rexp = 6.7 ± 
0.2 pc. The envelope has go = 15.3 ±0.7 Mq/pc^ and rexp ~ 
26.0 ± 0.3 pc. After 10 Gyr the values change to go = 750 ± 
15 Mq/pc"^ and rexp = 7.1 ± 0.1 pc for the core and go — 
6.6 ± 0.5 Mo/pc^ and rexp ~ 27.8 ± 0.8 pc for the envelope. 
Fig. ^ shows the data points for f = 1, 5 and 10 Gyr and the 
fitted functions for t = 1 and 10 Gyr. Plotted as vertical lines 
are the approximate tidal radii (Eq. The tidal extension 
of the MO beyond the tidal radius is clearly evident. 

The surface density is well fitted by a single de- 
Vaucouleur profile, 



E(r) 



Ev ■ exp -7.67 



1/4 



(13) 



with Ev = 1500 Mo/pc^ and = 18 pc (fitted at f = 
1 Gyr) (Fig. ||). Fitting the profile at t = 10 Gyr with two 
exponentials, results in a central surface density of 10600 ± 
300 Mq/pc'^ and exponential scale lengths of 10.8 ±0.4 and 
35.7 ± 0.5 pc, respectively. 

Measuring the line-of-sight velocity dispersion shows 
that the central velocity dispersion is about 13 km/s at the 
beginning and drops down to about 11 km/s after 10 Gyr. 
The velocity dispersion profile can be fitted by two expo- 
nentials (Fig. ^). 

The central density, surface density and velocity dis- 
persion decrease only very slowly and linearly with time 
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Figure 9. Evolution and destruction of the low-mass merger object R06. Shown are contour-plots with 20 pc resolution. Snapshots are 
at t = 200, 500, 700, 1000 and 1500 Myr. Highest surface-brightness (adopting AI/L = 1) of the MO is indicated above each panel. 



(Fig. g). The merger object is thus massive and stable 
enough to survive for a Hubble time despite the very large 
orbital eccentricity. 



3.3 Low-concentration super-cluster on an 
eccentric orbit 

If the initial super-cluster has a sufficiently small a and is 
subject to a strong tidal field, the merging efficiency can be 
low. The resulting merger object has a smaller mass and 
is less concentrated than model EOl above. Merger objects 
like this are easily disrupted at their first perigalacticon pas- 
sage forming an unbound stream of stars and surviving star 
clusters moving along almost identical orbits. However, a 
density enhancement survives for a substantial time, being 
composed of a core of particles with similar phase-space co- 
ordinates that may perhaps be identified as a dSph satellite 
(Kroupa 1997; Klessen & Kroupa 1998). 

In a calculation of a model with a Plummer radius of the 
super-cluster of 300 pc orbiting on an eccentric orbit with 
e = 0.33 only a small merger object R06, made out of 4 
star clusters, forms. It has a mass of about 0.35 • lO'' M© 
and gets dissolved after the first perigalacticon passage 
(Fig. m . The 'central' density drops from 100 Mq /pc^ down 
to 0.01 Mq /pc^ . The stars of this density enhancement cover 
an area a few kpc in diameter. Despite the fact that the stars 



do not form a bound entity, the density enhancement is still 
detectable after a few Gyrs. 



3.4 Evolution in the Kormendy diagram 

The models discussed in Sections 3.1-3^ are placed in the 
Kormendy diagram for a comparison with known stellar sys- 
tems (Fig. 0). 



The models on circular orbits (Section 3.1) overlap the 
Kr98 prediction shown as the six-pointed open star assum- 
ing M/Lb = 3. These models do not evolve noticeably over 
a Hubble time, apart from slight fading due to stellar evo- 
lution. 

The evolution of the dense super-cluster on a highly ec- 
centric orbit (Section 3.2) is shown for ages between 600 Myr 



and 10 Gyr as the solid curve in the thick box labelled EOl. 
Again, the evolution due to tidal harassment is not signifi- 
cant. The position of the equivalent 5 Myr old stellar cluster 
with Rpi = 50 pc and M = 1.23 x 10^ Mq is shown as the 
upper thick dotted open square labelled with the mass-to- 
light ratio M/Lb = 0.0144 (fig. 7 in Smith & Gallagher 
2001) and assuming B — V = 0. Fading without morpho- 
logical evolution (Eq. |l]) to M/Lb = 3 places such a cluster 
at the end of the thick long-dashed arrow in the vicinity of 
the Kr98 prediction, while the actual evolution (photomet- 
ric and dynamical) places it in the open box which lies in 
the vicinity of the fading vector of knot S (thick cross). 
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Figure 10. Known stellar systems are shown as faint open symbols, apart from highlighted knot S (Fig. |^). Models EOl and R06 
are shown at an age of about 5 Myr as the thick dotted squares assuming M/Lb = 0.014. Fading without morphological evolution to 
M/Lb = 3 would place these at the tips of the thick long-dashed arrows. The actual dynamical evolution is plotted as solid curves 
assuming the merger objects have M/L-q = 3. 



The evolution of the low-density super-cluster on an 
eccentric orbit (Section 3.3) is shown for ages spanning from 
200 Myr to 1.5 Gyr as the thick solid Une labelled R06. 
By 1.5 Gyr a remnant object stabilises in the vicinity of 
the known dSph satellites, and sports a population of 16 
co-streaming globular clusters. An equivalent stellar system 
that has the same Plummer density profile and mass as the 
initial super-cluster R06 appears at an age of about 5 Myr 
at the position of the lower thick dotted box labelled as 
M/Lb = 0.014. It would resemble a low-mass dirr, since 
the distribution of brilliant young star clusters would appear 
patchy in such a system. Hunter, Hunsberger & Roye (2000) 
discuss observable signatures of tidal-dwarf galaxies which 
may be related to this class of object. Fading to M/Lb = 3 
without morphological evolution would place this object at 



the end of the thick long-dashed arrow, in the vicinity of 
bright dSph satellites. 

Finally, from Fig. |l] we note that knot 430 (age « 
10 Myr), which is less massive than S but has a similar 
central surface density, may evolve via fading slightly past 
knot 225 (age ^ 500 Myr) to M/Lb = 3 if it undergoes 
no morphological evolution (eqn hi). Since the Milky- Way 
globular cluster ui Cen also lies in the vicinity of the final 
position of faded knot 430, we may assume that the precur- 
sor of Lo Cen may have looked similar to knot 430, or an 
object even more massive lying closer to knot S if dynami- 
cal evolution is taken account of (Fig. ^o|). uj Cen may thus 
have been born as a super-cluster containing hundreds of 
star clusters. Such a scenario may also be consistent with 
the unusual metallicity distribution in lu Cen (e.g. Smith et 
al. 2000), which is being taken as evidence that this unusual 
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cluster may in fact be a compact dwarf galaxy. Such an 
object encompassing about one hundred pc most probably 
forms over 5-10 Myr (Elmegreen et al. 2000) so that self- 
enrichment through supernovae is very likely. This scenario 
merits further study. 
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